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1. Nuclear shadowing or the depletion at low Bjorken x of the nuclear Deep Inelas- 
tic (DI) structure function, F%, with respect to the nucleon one, F% , has been observed 
in a number of experiments ( and references therein). Assuming the universality of 
parton distributions in nuclei, one expects nuclear shadowing to be present in other high 
energy processes as well, such as Drell-Yan pair, J/ip and T production in lepton-nucleus, 
hadron-nucleus and nucleus-nucleus collisions ( fU and references therein). In particular 
nuclear shadowing might be concurring with the other mechanisms, among which is the 
quark-gluon plasma formation, that result in a depletion of the observed cross sections for 
these processes. Moreover, as the very low x regime has become accessible at HERA, new 
phenomenological studies of high density QCD at the saturation scale predicted in are 
now possible ||. A quantitative understanding of both the x and Q 2 dependences of the 
nuclear parton distributions at low x therefore constitutes a practical and necessary step 
both for interpreting the outcome of future experiments at RHIC and at the LHC and for 
investigating the onset of parton saturation. 

Recent calculations rely on non-perturbative models for the nuclear parton distributions 
at a given (low) scale, Q 2 , combined with DGLAP || perturbative evolution. They are all 
therefore affected by the uncertainty in the initial parton distributions and, in particular, 
in the gluon distribution which governs evolution at low x, and which is poorly known 
experimentally. A strong effect is seen by changing the value of Q 2 itself which can plausibly 
vary within the range, Q 2 Q = 0.8 — fewGeV 2 , leading to sensibly different values for the 
shadowing of both the structure function and the gluon distribution at large Q 2 . 

Now, as Q 2 increases the (low x) gluon distribution in a proton should tend to a universal 
asymptotic value, corresponding to the Double Leading Logarithmic Approximation (DLLA) 
result of Ref. 0. Numerically, this value is attained at Q 2 ^ 10 3 GeV 2 provided the initial 
distributions grow much slower than « x~ 1 ^ 2 . Because of the coupling between the singlet 
quark and gluon distributions' evolution, a similar behavior is predicted for the structure 
function, F2 |J7|. It is therefore natural to address the question of whether the differences 
in the initial non perturbative nuclear shadowing will decrease with growing Q 2 and give 
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rise to a universal asymptotic curve which is within the range of planned experiments. Our 
first result is a negative one: we will show that because of the form that DGLAP evolution 
takes in a nucleus the influence of initial conditions is carried on to the largest attainable 
Q 2 values. 

We then examine carefully the asymptotic behavior of the shadowing ratios Rq = Ga/Gn 
and Rf = F^/FJ? , in the DLLA (notations are: Gtv(a) and F^^ for the gluon distributions 
and the structure function in a nucleon, (N), and in an isoscalar nucleus, (A), respectively). 
Our goal is to ascertain whether it is possible to distinguish among the different approxima- 
tions the perturbation series takes in a nucleus and at very low x. 

The proton structure function data analyzed recently at HERA |8||| have been shown 
|8j to lie in the asymptotic region (Q 2 3> Q 2 = 1 GeV 2 and i<i =0.1, and to evolve 
according to DLLA [[UJ. The key test is to prove that the data obey Double Asymtptotic 
Scaling (DAS) in the variables p = 7 ((F - Y o )/0 1/2 and a = ^((Y - F o )0 1/2 , 7 = 6/(33- 
2N f y/ 2 , Y = Inl/ac, f = ^ 2 \n(\nQ 2 /A 2 QCD / \nQ 2 JA 2 QCD ). Violations from DAS (other 



than due to the fact that the data lie in a pre-asymptotic region [TIJ) would signal either 
the onset of contributions beyond standard pQCD evolution ||, including the beginning of 
parton saturation 

The approach to asymptotia in a nucleus can be first analyzed by assuming that evolution 
proceeds through DLLA equations as well. We have found that in this case the ratios Rf 
and Rq display exact scaling in a, thus becoming a function of p only. As in the proton case, 
this is a model independent result in that we obtain a scaling form, independent from the 
initial conditions. We then use this result as a basis for addressing the next question i.e. the 
detection of violations of DAS scaling in nuclei, which could possibly originate at different 
values of x = x^ and Q 2 = Q 2 A , than in the proton. In particular Unitarity Shadowing 
Corrections (USC) |TT] , p!2| , pl| are expected to affect evolution at x^ > x Q because of the 
increase of the tranverse gluon density in a nucleus due to the overlapping of nucleons in the 
longitudinal direction. On a more speculative basis one might also expect the transition to 
the ln(l/x) resummation to appear in a different regime or, in the most "exotic" scenario, 



that medium modifications of the anomalous dimensions could be observed. In our approach 
such questions can be addressed systematically as they introduce specific scaling violations 
from the DLLA result, appearing as different a dependences in the ratios Rq and Rf- 

Our main observation is therefore that although it is technically predictable that for a 
proton target DGLAP evolution and DLLA should break down at very low values of x and 
sufficiently large Q 2 and give way to ln(l/x) summation and to USC, it is still a major 
task to be able to pinpoint where and if the transition from the different regimes is going 
to take place in the kinematical regimes currently under exploration. Our goal is to obtain 
some new insight by using nuclear targets where the asymptotic regime can in principle be 
reached at larger x. As a by-product we obtain quantitative predictions in the asymptotic 
kinematic regime which should be attainable at RHIC and at the LHC. 

2. We first summarize results for ordinary DGLAP evolution applied to the nuclear 
ratios at low x, assuming that the proton and the nuclear distributions evolve similarly. As 
it is well known evolution is driven by the gluon distribution which dominates over the sea 
quarks one and one can predict the behavior of the shadowing ratios, Rq and Rp with Q 2 : 

9R G f 1 r, fx ,„ 2 \\ G N {y,Q 2 ) r 2 , r*2<\ d V 



[ Pool-, * S (Q 2 ) Pf^ [Rg(v, Q 2 ) ~ Rg(x, Q 



d\nQ 2 Jx \y' J G N (x,Q 2 

__ OG N /dQ 2 {dG A {x,Q 2 )/d\nQ 



y 



-R G (x } Q 2 )) } (1) 



G N \dG N (x,Q 2 )/d\nQ 

dR F f 1 n> ( X /^2n\ G N (y,Q 2 ) r 2 ^2<\ d V 



where we have disregarded the sea quarks distribution on the r.h.s. of the coupled DGLAP 
evolution equations; P q c and Pqg are the splitting functions evaluated at NLO; and we 
used the approximation ~ 5/18Sjv(a)> with S = J2i Q.i{ x i Q 2 ) + <ji( x ,Q 2 )- F° r ease of 

presentation we will use the following notation: G' N ^ = dGjsf^(x, Q 2 )/d\nQ 2 . Eqs.(l) and 
(2) show that the Q 2 dependence of the ratios Rq and Rf is determined by a subtle balance 
involving both the parton distributions and the ratios themselves |T^]. Based on Eqs.(l) 
and (2), and defining Rg{x, Q 2 ) = R° G and Rf{x, Q 2 ) = R° F for the initial distributions, one 
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can make the following predictions for the Q 2 dependence of Rg and Rf'- i) Rg grows with 
Q 2 . In fact, if as predicted by non-perturbative shadowing models Rg is a growing function 
of x, then it also grows with Q 2 , the r.h.s. of Eq.(l) being positive (y > x); ii) if R° F < R G , 
then Rp grows with Q 2 ; if R° F > R G , then Rp initially decreases with Q 2 until it reaches 
the value of R G and it subsequently starts increasing along with Rg- 

Note that from the behavior of Rf and Rg obtained from the straighforward application 
of DGLAP equations, one cannot predict the approach of e.g. Rg to a universal limiting 
curve at large Q 2 . As a matter of fact, although the form of Eq.(l) might seem suggestive of 



a fixed point behavior |T5|], this is a priori not the case, since the quantity Ga'/Gn 1 depends 
on Q 2 . The rate of change with Q 2 is instead governed both by: (a) by the ratio, G' n /Gn] 
(b) by the difference A G = G' A /G' N — Rg(x,Q 2 ) at Q 2 = Q 2 . Current parametrizations 
|T6| , p!7|| feature an ultra-soft behavior of the gluon distribution at Q 2 < 1 GeV 2 , i.e. G^ — > 
as x — > 0, thus causing (because of (a)) a rapid evolution which strongly reduces the 
shadowing in both Rg and Rf, between Q 2 Q and Q 2 ~ 1 — 2 GeV 2 . If on the contrary 
one assumes Q 2 Q ranging from 2 to 5 GeV 2 , where harder low x initial gluons are expected, 
then the evolution is slower (GV is larger) and the nuclear ratio is basically unchanged 
at Q 2 ~ 10 — 100 GeV 2 . A further model dependence follows from the usage of different 
non-perturbative shadowing mechanisms. We examine two in particular: the Aligned Jet 
Model (AJM) (see JL8| and references therein), and Initial State Recombination (ISR) ||19|| . 



Both models explain qualitatively the initial onset of shadowing. Accurate quantitative 



calculations have been performed using the AJM in p0|j21[| . 

Results are summarized in Fig.[l], where we show the Q 2 dependence of the ratios 
R F (x,Q 2 ) and Rg(x,Q 2 ) in 40 Ca at fixed x = 10" 4 (Fig. 1(a)) and x = 10" 2 (Fig. 1(b)), 
for both the AJM (full lines) and the ISR model (short dashes). In order to show the de- 
pendence on the initial scale Qq, results are presented for both models by taking, Ql = 0.8 
GeV 2 , and Qq = 5 GeV 2 (the latter can be easily distinguished in the graph by observing 
the shift in the starting point of the curves). Moreover, as the main purpose of the figure is 
to illustrate the main features of both the nuclear models and the initial parton distributions 



that will lead to our description on the asymptotic behavior, we have not sought for Rp the 
best agreement with the data. Details on this part are going to be given elsewhere ||22|| . A 



remarkable feature is the sensitivity of the ISR model to this initial scale Qq. It is only for 
Qo ~ 1 GeV 2 that a sizeable gluon shadowing is obtained due to the fact that the amount 
of initial shadowing is proportional to the square of the initial gluon distributions and to 
a s{Qo) / Qo El]- R- esu hs obtaind with the AJM vary less dramatically with the initial scale, 
Qq, which in this case enters just the (qq) — nucleon (or (gg) — nucleon) cross sections |18 



Nonetheless, the initial difference is carried on to Q 2 as large as 10 3 GeV 2 . Moreover, in both 
cases we can observe the Q 2 behaviour outlined before: a rapid suppression of shadowing in 
the range of Q 2 up to 2 GeV 2 and a subsequent softer evolution. The comparison between 
the Q 2 behaviour for different fixed x values (Fig. 1(a) and Fig. 1(b)) shows that evolution 
is slower for smaller x. This can be technically understood by noting that, in the DLLA 
limit the logarithmic derivative of Rq in Eq. ([!]) vanishes. This is why initial discrepancies 
between models are more likely to persist at smaller x. 

In summary, the asymptotic behavior of Rq depends on the initial conditions up to the 
largest values of Q 2 attainable at low x. This feature is in common with the proton gluon 
distributions themselves as shown e.g. in [0. As it is well known, neither the data nor 
theoretical arguments can help us defining the optimal values of Rq(x, Q 2 ) and Rp(x, Q 2 ). 

The situation that we have described calls for some redefinition of the approach to 
asymptotia in deep inelastic scattering from nuclei. In the next Section we illustrate how 
different behaviors of the data could be revealed by extending the double asymptotic scaling 



analysis of Ref. |T0| to nuclei. 

3. We now examine the nuclear DI structure function and gluon distribution in the 
asymptotic regime defined by x — > 0, Q 2 —>■ oo. Our goal is to explore scaling relations 
in nuclei in order to be able to compare theory with data in a model independent way. 
The derivation of the equations of the DLLA in a nucleus parallels the one for the proton, 
namely one first writes the DGLAP evolution equation for the gluon distribution in the limit 
n — > 1, n being the variable in moments space (we have omitted the subscripts N(A) unless 
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necessary) : 

dg(n,Q 2 ) a s (Q 2 ) , 

-TGG( n )9{n,Q ), (3) 



<91nQ 2 2tt 

Jgg(n) ~ 2Ca/ (n — 1) + k being the anomalous dimension in the limit n ^ 1 (n = — 11/6 — 
rif/3CA is the next-order or subleading contribution in this limit). Solutions in (x, Q 2 ) are 
found by evaluating the anti-Mellin transform, 

G(x, Q 2 ) = ^~J c dn g{n, Q 2 ) exp [Y(n - 1)] , (4) 

with the saddle point method |§. In Eq.(|), g(n,Q 2 ) = g(n, Ql) exp [£/(n - 1) + £«/2]; 

Ql) = Jq dx x^ n ~^ g(x, Ql), g(x, Ql) being the initial gluon distribution, and G(x, Q 2 ) = 
xg(x,Q 2 ); Y = ln(l/x), and £ = 7 2 ln(ln(Q 2 /A 2 )/(hiQ 2 /A 2 )), ^ = c A /(irb), b = (33 - 
2Nf)/12ir. We rewrite the integrand in Eq.(f|) as: 

g(n,Q 2 ) =g(n,Q 2 )ex P ^/2ex V [Yf 1 (n)+^f 2 (n)], (5) 

where fi{n) — n — 1, /2(^) — (n — l)" 1 and K and £ are both similarly large. If one takes 
"soft" initial conditions such as, G(x, Ql) ~ A N x~ x , A < 0, then Q 2 ,) ~ A N /{n — (A + l)) 
has a pole to the left of the saddle point which is found by setting d[Y fi(n)+£f2(ri)]/dn = 0. 
The integral in Eq.(f|) is then approximated by: 

^ nr ( 9(n ,Q 2 ) X ' " 



^ gKQ >xp(W2) ^^-tiM^ (6) 

l*7iK) +^K)| 



where n = 1 + y^/AY, AY = Y — Y a , Y a = ln(l/x D ), x Q « 0.1, is the saddle point and 
/f( 2) K) = d 2 /i(2)(n)/dn 2 [ =no . 

By introducing the variables, p = j((Y - Y )/^) l/2 and a = 7~ 1 ((V - F G )£) 1/2 , Q one 
has: 

n = n„(p) = 1 +j/p, (7) 

and, 
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G = G-V ,) = Mh ) { 1) gfeM , (8) 

where fc{p/l) is a smooth function describing the initial conditions. 

DAS is the prediction that, in the hypothesis of soft initial conditions, and in the asymp- 
totic limit defined by o — > oo and p « 0(1), \a(G DAS / fa(p, c)) becomes a linear function 
of a, independent of the value of p, with slope fixed by the known constant, 7. 

A similar behavior is found for Ff JTOf and can therefore be compared to the available 
data |8|||]. The structure function's asymptotic behavior is in fact obtained by solving the 
equation [} 

dF$(x,Q 2 ) 5«» . ^ 2 



<91nQ 2 18 vr 

yielding: 



G(x,Q 2 ), (9) 



^(p^) = /s(p,a)exp(2 7( T), (10) 

where /s was derived by using the LO expression for a s and, as for the gluons, it depends 
on the initial conditions. The scaling of F% can be seen from Fig.^ (top-right). From Fig.|2| 
one can also see that the data from NMC (triangles) p^j , corresponding to larger x with 
respect to the 1995 HERA ones (open squares) ||, as well as some of the more recent HERA 
data with very low x and Q 2 (open dots) 0, violate scaling (for a better reading compare 
with top-left). These scaling violations have been interpreted as due to the fact that the 
kinematics is not yet asymptotic, as it can be easily seen from the fact that the data lie well 
below <t ~ 1. 

Deviations from DAS can be predicted also in the case of hard initial conditions i.e. 
when g(n ,Ql) w Aisf/(n — (A + 1)) with A ^ 0.2. Intuitively this corresponds to taking 
the limit Y ^> £ in Eq.(f|), thus defining the new saddle point: n Q = (1 + A) + 1/AY. The 
corresponding gluon distribution is then 



1 We refer here to the singlet part of the structure function 
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G(x, Q 2 ) = f h G exp [X(Y - Y ) + e/A + £k/2] 
= /£exp [Aap + ( 7 2 /A + k/2)£ 



(11) 



with flj(p,a) = (see also [0]). This behavior supports the presence of USC 



appearing as a non-linear term in the evolution equation [IT], which has the effect of damping 
the steep rise of the gluon distribution at small x. An alternative possible explanation of 
DAS violations in the recent HERA data || which accounts also for the peculiar stooping 
of the logarithmic slope of the proton structure function, dF 2 /dlnQ 2 , at low x and Q 2 , is 
that indeed USC |24j] need to be taken into account. 

We now study these two different scenarios for the asymptotic behavior in nuclei, where it 
is well known that some aspects of perturbative evolution such as USC, are expected (simply 
based on geometrical arguments) to arise at larger values of x. In Fig.0 we present the world 
low x data on the nuclear ratios, Eq.(|l|) as a function of the DAS variable, a (bottom-right), 
and we compare both data and their kinematics (bottom-left) with the proton ones (top). 
We use in both cases the following values of Q 2 , Aqcd, P and a: Q 2 = l GeV 2 , Aq C d = 185 
MeV (LO). f\ From the figure it appears that the existing nuclear data are scarse and do 
not presently support a DAS type analysis and, moreover, they seem to lie mainly in a 
pre-asymptotic region. Experiments at RHIC and LHC are however expected to be able to 
cover the asymptotic region. Although experimental extractions of the logarithmic slopes in 
nuclei have been performed in Jl], very little can be concluded from these data as well ||22|| . 

We evaluate the ratios Rq and Rp, Eq.(|l]) in DAS, i.e. within the hypotheses: i) the 
quark and gluon distributions are initially shadowed due to some non-perturbative mech- 
anism; ii) the pQCD evolution mechanism is not affected by the nuclear medium; in) the 
initial distributions are soft. The DLLA predictions are: 



R G (x,Q 2 ) 



gA(n ,Q 2 
9N(n ,Q 2 



3/2 



[ 9l(n ,Q 2 



1/2 



(12) 



2 In the NLO analysis performed in [1C] it was found that: Q 2 = 
(NLO), pss 1 and a > 1.2. 
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1.8 GeV 2 , Aqcd = 200 MeV 



where n Q is the saddle point defined by Eq.(^) for both the proton and the nucleus We 
rewrite Eq. (|1^) as a function of the DAS variables by using Eq. (|8|) : 

r>DAS . 

Re -jsm (13) 

i.e. the exponential terms appearing in G DAS (p,a), Eq.(|3|), are the same in a nucleus and 
in a single nucleon respectively, thus canceling the a dependence in the ratio Rg'- Rg is 
predicted to scale exactly in a to a smooth function of p that is determined entirely by the 
(soft) initial conditions. 

The onset of a different evolution mechanism in the nucleus will appear as a a- scaling 
violation modifying the exponential behavior of Eqs.fl3D-(|l"2"D with respect to the single nu- 
cleon case. Since the low x behavior of recent HERA data seem to show evidence for rather 
large screening corrections, and at the same time they do not rule out the hard pomeron 
contribution, we consider the effect of: (A) USC combined with soft initial conditions; (B) 
USC combined with hard initial conditions; (C) Hard initial conditions in both nucleon and 
nucleus, no USC. 

The effect of USC is taken into account through a "damping factor", Dg(x,Q 2 ) < 1, 
(Ref. |24| and references therein), evaluated using Mueller- Glauber's eikonal approximation 
|27fl . We extended the calculation to nuclei by assuming that two gluons inside a nucleus are 
correlated by a larger confinement radius, Ra ~ tqA 1 ^ 3 , than in a nucleon, corresponding 
to a smeared impact parameter space two-gluon form factor (details of this calculation will 
be included in [|^| ). As a result the effect USC is enhanced in a nucleus with respect to a 
nucleon target, due to the larger transverse gluon density at similar values of x. 

The asymptotic gluon distribution function is written in terms of the damping factor as: 

Gf% A) (p,v) = DS {A) x&> A3 M. (14) 
By using Eq.(|l4|) We can now calculate the ratio Rg, for the cases listed above. We obtain: 



R^\p,a) = R D G AS {p) 



x 



exp 


2j(a 


- ca) - 


7 2 K cr 
2 p_ 


exp 


2 7 (<7 


- ojv) - 


2 ~P_ 
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R% AS (p) x exp [ {c? a - <tn)\ ■ 



(15a) 



R ( G B \p,a) 



h,A 
G 



exp 



JG 



X 



exp 



AatCTP + 



7 £ 



Oat 



= C exp 



(A A - Ajv)crp + 7 2 ( 



■* (p, a) = C exp 



A, 
1 



Aat p 
1 ,cr 



x exp [-{a a - vn)] ■ 



(A A - Ajv)ap + 7 (- — )- 



(15b) 
(15c) 



Aa A ^ p 

Here Rq AS is the same as in Eg. (pT3|) ; (Tn(a)(p, c) = l/27 m (l/.DGr); an d the shadowing 
for the initial hard distributions has been parametrized as Rq = f G A j '/q N = Cx a , with 
C ~ 1.3 and a ~ 0.08 — 0.1. Moreover, we have chosen Q 2 Q = 1 GeV 2 for both soft and 
hard initial conditions and A at = 0.35. Our scaling result, Eq . (|T3|) , and the scaling violating 
ones, Eqs.(15a)-(15c), are shown in Fig. ^| for two different values of p: p = 1.8 well inside 
the asymptotic region shown in Fig|J, and p = 3.4 corresponding to very low x and almost 
fixed Q 2 where we expect standard DGLAP to break down. 

A few comments are in order. Starting from soft initial conditions one obtains either the 
cr-scaling curves (full lines) or the a-scaling violation, Eq.(15a), induced by USC (dashed 
lines). These corrections are driven by the damping factor given for each value of p by the 
dashed curves below. The decreasing trend with growing a is larger at large p because of the 
correspondingly decreasing values of x at similar values of Q 2 . On the other side, hard initial 
conditions, Eq.flTTp and dot-dashed curves in Fig.3, show sensible deviations from a-scaling 
at large p where the ln(l/a;) term dominates over the InQ 2 one. This effect is enhanced 
by USC (dotted curves). An interesting observation is the change in slope of the scaling 
violations both when passing from soft to hard initial conditions at low p, and when passing 
from low p to large p. These and similar other regularities could be studied systematically 
both for the gluon and the structure function ratios once a much larger set of data will be 
available. 

Most importantly, with the approach proposed here we eliminate the ambiguities in the 
determination of the value of gluons and quarks nuclear shadowing illustrated in Fig.l, 
in that we identify scaling relations that must be verified independently from the initial 
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non-perturbative nuclear shadowing. 

In conclusion, by applying DAS to nuclei we have shown model independent predictions, 
i.e. scaling in the variable p for the ratios of the nuclear gluon distributions to the free 
nucleon ones in the asymptotic region. Similar relations hold for the nuclear structure 
functions. We have considered a few possible sources of scaling violations due to the onset of 
USC and to the domination of a hard pomeron. As a result, by studying the A-dependence of 
shadowing we find some new constraints on perturbative evolution at low x. Our calculations 
are relevant for the regime accessible at future experiments at RHIC, LHC and at the eA 
project at DESY [ f2~8|1 . 
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FIG. 1. DGLAP evolution in 40 Ca. Experimental data from NMC |25| (triangles) and E665 
(stars); theoretical curves are explained in the text. 
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FIG. 2. The asymptotic region in the proton (top) and in a nucleus (bottom). The figures on 
the left show the kinematical range in Y = 1/x and £ covered by current experiments: NMC |23| 
|| (triangles), HI §a (squares) and §b,c (full dots), ZEUS §§ (open dots), E665 @ (stars). 
The asymptotic region is delimited by the values: 0.7 < p < 3 (full lines) and 1 < a < 2.1 (dashed 
lines). The figures on the right show the onset of DAS in the proton (top) and the shadowing ratio, 
Rf, vs. a in a nucleus (bottom). While DAS is achieved in the asymptotic region covered mostly 
by the HI || data, a regular pattern is yet to be seen in the nuclear data which are clearly lying 
outside the asymptotic region (left bottom). 
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FIG. 3. Top: Ratio Rq = Ga/Gn, vs. DAS variable a at two different values of p: p = 1.8 (left) 
and p = 3.4 (right). The theoretical curves show both cr-scaling (full lines) and the cr-scaling viola- 
tions predicted in this paper: dashes, Eq.(15a), dots, Eq.(15b), and dot-dashes, Eq.(15c). Bottom: 
Ratio of nuclear to nucleon damping factors (see Eq. (|T4|) ) , as a function of x for the same values 
of p and a used in the shadowing calculations above. The Q 2 ranges are: 2.3 < Q 2 < 2 x 10 3 GeV 2 
(left), 1.5 < Q 2 < 40. < GeV 2 (right). We show results using the AJM of Ref. (T^j as our 
non-perturbative input; analogous scaling relations are found by using other initial shadowing 
models. 
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